Spin and charge dynamics in [TbPc 2 ]° and [DyPc 2 ]° single molecule magnets 
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Magnetization, AC susceptibility and fiSK measurements have been performed in neutral ph- 
thalocyaninato lanthanide ([LnPc2] ) single molecule magnets in order to determine the low-energy 
levels structure and to compare the low-frequency spin excitations probed by means of macroscopic 
techniques, such as AC susceptibility, with the ones explored by means of techniques of microscopic 
character, such as fiSK. Both techniques show a high temperature thermally activated regime for 
the spin dynamics and a low temperature tunneling one. While in the activated regime the cor- 
relation times for the spin fluctuations estimated by AC susceptibility and /iSR basically agree, 
clear discrepancies are found in the tunneling regime. In particular, /xSR probes a faster dynamics 
with respect to AC susceptibility. It is argued that the tunneling dynamics probed by /xSR involves 
fluctuations which do not yield a net change in the macroscopic magnetization probed by AC sus- 
ceptibiliy. Finally resistivity measurements in [TbPc2] crystals show a high temperature nearly 
metallic behaviour and a low temperature activated behaviour. 

PACS numbers: 75.50.Xx, 76.75.-K, 76.60.Es 



I. INTRODUCTION 

The trend towards ever-smaller electronic devices is 
driving electronics to its ultimate molecular-scale limit, 
which allows for the exploitation of quantum effects. 
Single Molecule Magnets (SMM) are among the most 
promising materials to be used in molecular spintronic 
devices^ or as logic units in quantum computers^, since 
they combine the classical macroscale properties of bulk 
magnetic materials with the advantages of nanoscale en- 
tities, such as quantum coherence. It has already been 
shown theoretically^ that molecular magnets can be used 
to build efficient memory devices 3 - and, in particular, one 
single molecule can serve as a storage unit of a dynamic 
random access memory device^. In addition, SMM have 
been recognized as suitable materials which can be used 
as contrast agents for medical diagnostics. In fact, these 
molecules comprise the chemical and magnetic structure 
of the gadolinium chelates with the superparamagnetism 
of the SPIOs (Superparamagnetic Iron Oxides), materials 
which represent some of the most common MRI contrast 
agents in current use£. 

LnPc2-based compounds (Pc=C32Hi6Ng is phthalo- 
cyanine, Ln a lanthanide ion) represent a new class of 
SMM characterized by a doubly-degenerate ground state 
and large magnetic anisotropy which yields a remarkably 
large barrier to spin reversal. In particular, [TbPc2]~ 
and [DyPc2]~ complexes are the first mononuclear sys- 
tems exhibiting a very slow relaxation of magnetization 7 
and consequently a characteristic correlation time (t c ) 
for the spin fluctuations reaching several /is at liquid ni- 
trogen temperature. At variance with transition metal 
molecular magnets, in the LnPc2 molecules the electronic 
levels of the ground-state multiplet, with angular mo- 
mentum J, are mainly split by the strong anisotropy of 



the crystal field (CF) at the Ln 3+ ion. Accordingly the 
separation between the double degenerate ground state 
and the first excited levels can reach several hundreds 
of kelvin, thus yielding unprecedently large r c at cryo- 
genic temperatures. Moreover, it is noticed that the 
magnetic properties of these lanthanide complexes are 
strongly influenced by the chemical and structural mod- 
ification of the ligands through which the CF potential 
can be varied^. The CF configuration can be modified 
from the very dilute limit in [LnPc 2 ]~[TBA]+N[TBA]Br, 
with N+lS> 1 the number tetrabutylammonia units, to 
the concentrated limit consisting solely of neutral double- 
decker molecules [LnPc 2 ]°. Recently it has been shown 
that one or two-electron oxidation of [LnPc2]~ yields a 
compression of the cage around Ln 3+ ion and an increase 
of the CF splitting^*! 1 .. 

By means of nuclear magnetic resonance (NMR) and 
muon spin relaxation (/iSR), the correlation time for the 
spin fluctuations in the neutral compounds [TbPc2]° and 
[DyPc2]° was found to be close to 0.1 ms at 50 K 1 ^, 
about two orders of magnitude larger than the one previ- 
ously observed in the non-oxidized lanthanide based sin- 
gle molecule magnets. In [TbPc 2 ]° two different regimes 
for the spin fluctuations have been evidenced: a high tem- 
perature thermally activated regime involving spin fluc- 
tuations, driven by spin-phonon coupling, across a bar- 
rier A ~ 880 K separating the |J = 6,m = ±6) ground 
states and the \J = 6,m = ±5) first excited states, 
and a low temperature regime involving quantum fluctu- 
ations within the twofold degenerate ground-state. For 
[DyPc2]° a high temperature thermally activated regime 
is also found, however it cannot be explained in terms of 
a single spin-phonon coupling constant. 

Here we report an analysis of the static uniform spin 
susceptibility (xs) f° r the [DyPc2]° complex, from which 
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a determination of the low energy level structure of the 
Dy 3+ ion spin multiplct J=15/2 is made. On the basis of 
these results for the CF level splitting it is possible to ana- 
lyze the temperature (T) dependence of muon spin-lattice 
relaxation rate in the presence of an applied held, de- 
rived from previous /iSR measurements performed at the 
ISIS facility 12 and from new measurements performed at 
the Paul Scherrer Institute (PSI). Accordingly, the spin- 
phonon couplings driving the transitions among the CF 
levels have been derived and the characteristic correla- 
tion time for the spin fluctuations estimated. /iSR results 
are compared to the ones obtained by AC susceptibility 
measurements in the same samples. It is found that both 
AC susceptibility and /zSR probe a high temperature ac- 
tivated and a low temperature tunneling regime for the 
spin dynamics. However, while the values estimated for 
the correlation times in the high temperature activated 
regime are quite close, at low temperature the tunneling 
rate probed by /iSR is much faster than the one derived 
by AC susceptibility. Finally, transport measurements in 
[TbPc2]° crystals are presented. 
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FIG. 1: Temperature dependence of xT in [DyPc2]° for 
H = 1000 Gauss (open circles) and the best fit curve cal- 
culated according to Eq.l curve (solid line). The dashed 
line corresponds to the behaviour calculated on the basis of 
[DyPc 2 ] _ TBA + CF splittings derived by Ishikawa et al. [H 



II. TECHNICAL ASPECTS AND 
EXPERIMENTAL RESULTS 

All reagents were purchased from Across or Aldrich 
and used without further purihcation. The synthesis 
of both samples was carried out by using several mod- 
ifications of the protocol based on the published pro- 
cedure reported in Refill. The synthesis was accom- 
plished by templating reactions, starting from a mix- 
ture of the phthalonitrile precursor o-dicyanobenzene 
and the lanthanide acetylacetonate Ln(acac)3- nH^O, 



in the presence of a strong base (e.g. DBU, alkox- 
ides) and high-boiling solvents, such as pentanol or 
hexanol. A mixture of 1,2-dicyanobenzene (42 mmol), 
Ln(acac) 3 -4H 2 (Ln = Tb 3 +, Dy 3 +) (5 mmol), and 1,8- 
diazabicyclo[5,4,0]undec-7-ene (DBU) (21 mmol) in 50 
mL of 1-pentanol was rehuxed for 1.5 days. The solution 
was allowed to cool to room temperature and then acetic 
acid was added and the mixture was heated at 100 C 
for 0.5 hours. The precipitate was collected by filtration 
and washed with n-hexane and Et20. The crude purple 
product was redissolved in 800 Ml of CHCl 3 /MeOH (1/1) 
and the undissolved PcH2 was filtered off. Both forms, 
blue (anionic [LnPc2]~) and green (neutral [LnPc2]° ), 
were obtained simultaneously, as revealed by electronic 
absorption spectra^. In order to convert the unstabilized 
anionic form to the neutral one, the reaction mixture was 
presorbed on active (H2O-0%) basic alumina oxide. Pu- 
rihcation was carried out by column chromatography on 
basic alumina oxide (deactivated with 4.6 % H2O, level 
IV) with chlorophorm-methanol mixture (10:1) as eluent. 
In general, the yield was 30-35%. By means of additional 
radial chromatography on silica gel followed by recrys- 
tallization from chloroform-hexane mixture, analytically 
pure powder samples were achieved. 

Deep green crystals of the products were obtained by 
using slow diffusion of CH2CI2 into C2H2CI4 solution of 
the pristine [L11PC2] . After 2 weeks, deep green needle- 
like crystals were obtained (see Figj2]). The [L11PC2] 
molecules crystallized in the space group P2i2i2i (7- 
phase) as reported in ReffTol and were isomorphous to 
each other. 




FIG. 2: Photograph of a few [TbPc 2 ]° crystals grown accord- 
ing to the procedure reported above. The numbers in the 
ruler scale correspond to centimeters. 

DC Magnetization (M) and AC susceptibility mea- 
surements have been performed by means of an MPMS- 
XL7 Quantum Design SQUID magnetometer. The static 
uniform susceptibility was determined from the ratio 
Xs = M/H in an applied field H = 1000 G in the 2- 
300 K temperature range. In Fig. ^ the T dependence 
of xsT is shown. In-phase (x'-real) and out-of-phase (%"- 
imaginary) AC susceptibility data were collected at oscil- 
lating held frequencies between 10 and 1488 Hz, in the 10- 
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80 K temperature range and for different applied static 
fields. In Fig. ^ the temperature and frequency depen- 
dence of the AC susceptibility for [TbPc2]° is shown. A 
sharp drop in x'T is observed on cooling in correspon- 
dence to a peak in x" /xs, which progressively shifts to 
lower temperatures as the AC frequency is lowered. In 
[DyPc2]° analogous but slightly broader peaks are ob- 
served at lower temperature (Fig. @}, suggesting a faster 
dynamic (see Discussion). 
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FIG. 3: Temperature dependence of a) x'T and b) x" Ixs in 
[TbPc2]°, where x ', x" and xs are in-phase-AC, out-of-phase- 
AC and DC molar magnetic susceptibilities, respectively. The 
measurements were performed in 4 G oscillating magnetic 
field at 10, 300, 997 and 1488 Hz, in the presence of a 1000 G 
DC component. 

Longitudinal field /iSR measurements in the 50 mK- 
300 K temperature range were carried out either at ISIS 
pulsed muon facility on MuSR beam linei^, or at PSI 
muon facility on GPS (General Purpose Surface-Muon 
Instrument) and LTF (Low Temperature Facility) beam 
lines. In contrast to the ISIS measurements, the contin- 
uous beam muon production at PSI allows for a much 
more accurate determination of the short-time relax- 
ation. In Fig. (|5|) zero-field (ZF) depolarization curves 
at a few selected temperatures are shown for [DyPc2]°. 
A Kubo-Toyabe (KT) relaxation 1 ^ is clearly observed be- 
low T* ~ 60 K, indicating very slow fluctuations, that is, 
vj (7(n-\/ (A^Dm) *C 1, with v — 1/t c being the character- 
istic frequency of the spin fluctuations and (A/i^_) p is the 
mean squared amplitude of the field fluctuations at the 
muon site. 7^ is the muon gyromagnetic ratio. By fitting 
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FIG. 4: Temperature dependence of the normalized out-of- 
phase spin susceptibility x" Ixs in [DyPc2]°, for a 4 G oscil- 
lating magnetic field at various frequencies. A 9000 G static 
external field was applied. 
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FIG. 5: Time evolution of the muon polarization in zero-field 
in [DyPc2] sample normalized to its value for t — > at three 
selected temperatures. 



the low-temperature data with a KT function one can es- 
timate the static field distribution probed by the muons 
around y/{Ah^ = 490 G and y/{AK\)^ = 550 G for 
the Dy and Tb compounds, respectively. These values 
are in excellent agreement with the measured "repolar- 
ization" of the signals with applied longitudinal fields. In 
the slow fluctuations regime only the 1 /3 tail of the KT 
depolarization function is affected by slow fluctuations 
and it decays according to P^(t) = exp(—(2/3)vt) — 
exp(-Xt). On the other hand, for T > T* (around 60 
K for [DyPc 2 ]° and around 90 K for [TbPc 2 ] ) the fast 
fluctuations limit is attained. In this limit the decay of 
the polarization in zero or in a longitudinal field is given 
by P^(t) = expi^liAh^^t/v) 13 = exp(-XtY , where 
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the exponent j3 can differ from unity in case of a dis- 
tribution of muon sites, as it can be the case here. In 
the intermediate fluctuation limit v / {-^^^J (Ah 2 ^) ^ > 1 
(T~T*), A starts to decrease according to the so-called 
Abragam formic. 
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FIG. 6: T dependence of the muon ZF longitudinal relaxation 
rate (circles) in [TbPc2]° and [DyPc2]° (in the inset), together 
with 1000 G longitudinal field data from Ref. []3 (squares). 

By using the aforementioned expressions for the de- 
cay of the polarization in the fast and slow fluctuations 
regime it was possible to derive the T-dependence of A 
(Fig|6]) in zero-field and estimate the T-dependence of r c 
(see Fig. [7]), which was compared with the one previ- 
ously derived at ISIS for a 1000 G longitudinal field-^. 
The data for both materials were fit with a stretched 
exponential above T*, with a T-independent exponent 
P = 0.5. In Fig. [7] one notices that in the high T regime 
the t c derived in zero-field and in a 1000 Gauss longi- 
tudinal field are in satisfactory agreement and show an 
activated T-dependence of r c . On the other hand, below 
50 K, where tunneling processes become relevant, t c (T) 
is observed to flatten at values which differ by more than 
an order of magnitude upon increasing the field from zero 
to 1000 G. The increase of r c with the magnetic field in- 
tensity, also observed in the AC-susceptibility measure- 
ments (Discussion) , should be ascribed to the removal of 
the ground-state degeneracy by the magnetic field, which 
progressively inhibits the tunneling processes. 

Data taken below 1 K in LTF yield an essentially tem- 
perature independent depolarization rate for both com- 
plexes down to 50 mK, as expected for tunneling domi- 
nated depolarization. However, there exists a weak max- 
imum in A(T), near T = 0.2 K for the Dy sample and 
near T = 0.7 K for the Tb sample, which warrant further 
studies (Fig. ((5]) and see Discussion). 

Resistivity measurements on single crystal samples of 
[TbPc2]°, with current along the c-axis, have been per- 
formed in the 2.8 - 294 K temperature range in DC mode 
by using a four-terminal technique. Contacts were made 
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FIG. 7: T dependence of the correlation time for the spin 
fluctuations in [TbPc 2 ]° and [DyPc 2 ]° (in the inset, for T > 
T*) derived from A data reported in Fig. ©. 
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FIG. 8: T dependence of the muon relaxation rate in zero-field 
for [DyPc 2 ]° and [TbPc 2 ]°, at T < 1 K. 



by attaching 0.01 mm diameter gold wires to the samples 
with conductive paste and sub-^(A currents were used. 
The samples were left free-standing to minimize thermal 
stress and consequent cracking. The voltage has been 
measured with a voltmeter with an internal impedance 
larger than 10 Gf2. Three different crystals have been 
measured and they all showed an increasing resistance 
on cooling from 2.5 kohm to greater than 2 Mohms be- 
low 15 K. In Fig. (j^J) the temperature dependence of the 
resistivity for a crystal with a 2.72 • 10~ 3 area/lenght ra- 
tio is shown. The sample conducts at room temperature 
with a resistivity p ~ 6.4 f2-cm. Upon cooling, p presents 
an initial smooth decrease till about 220 K and then an 
activated growth, characterized by an activation energy 
AE ~ llmeV, down to about 25 K. 
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FIG. 9: Temperature dependence of resistivity p, measured 
for a [TbPc2]° crystal. Solid line represents the best fit of the 
activated trend in the 220 - 30 K region. 



III. DISCUSSION AND CONCLUSIONS 

First we shall discuss the temperature dependence of 
the static uniform susceptibility xs m [DyPc 2 ]°, in order 
to derive the CF levels structure. Since we are deal- 
ing with a powder, the molecules and the CF axis are 
randomly oriented with respect to the external field H 
| z. Thus the total magnetization M can be deduced 
by summing the contributions of the molecules with the 
anisotropy axis perpendicular and parallel to z, namely 
M(T) = |M S , B (T) + \M Z {T), with 



Mi{T) = N A 



V +15/2 (u k )e~ 



E 



+ 15/2 
fc=-15/2 ' 



E<< 



x,y,z 



(1) 



In Eq. (|TJ|, represents the k-th eigenvalue of the 
hamiltonian: 



energy splittings are Ai = 65 K, corresponding to the 
separation between the \m = ±13/2) ground states and 
the |m = ±9/2) first excited levels, A 2 = 47 K between 
the \m = ±9/2) levels and the \m = ±11/2) second ex- 
cited levels and A3 = 460 K between the |m = ±11/2) 
levels and the \m = ±15/2) third excited levels. These 
values differ from the ones deduced for [DyPc2]~ on the 
basis of crystal field calculations (Ai ~ 50 K, A2 oi 207 
K, A 3 ~ 126 K) 13 . As it can be seen in Fig. (JTJ, a good 
fit of the experimental data is found. 

From the CF splittings, combined with the /iSR relax- 
ation data (Fig. (TTU)) '), it is possible to derive information 
on the spin-phonon coupling driving the high T spin fluc- 
tuations in [DyPc2]°. Since the energy difference between 
the muon hyperfine levels and the m levels of Dy 3+ spin is 
large, the muon longitudinal relaxation rate A is driven 
by an indirect relaxation mechanism involving a muon 
spin flip without change in m. This is possible thanks to 
the tensorial nature of the hyperfine coupling constant, 
which allows the coupling of the transverse components 
of the hyperfine field h x _ y to J z . Thus, denoting with 
(r m ) the finite life-time of the crystal field levels induced 
by the spin-phonon scattering processes, A can be written 
in the formii 



+15/2 
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= -15/2 



-Em/T 



(3) 



J L'm 



E m being the eigenvalues of the CF levels and Z is the 
corresponding partition function. It is noted that the low 
magnetic field (1000 Gauss) applied during /xSR exper- 
iments yields a negligible correction to E m and, hence, 
its effect is negible for kgT HbH. The life-time for 
the m levels can be expressed in terms of the transition 
probabilities p m ,m±i between m and m ± 1 levels, which 
depend on the CF eigenvalues and on the spin-phonon 
coupling constant G^: 



1 



Pm,m — 1 H - Pm,rn-\-l i 



(4) 



H — Hcf + g^B J • H 



(2) 



where TicF is the crystal field hamiltonian and g/isJ ■ H 
is the Zeeman term, (/jf?) is the expectation value of 
the i-th component of the magnetic moment over the k- 
th eigenstate of the hamiltonian. In order to analyze 
XsT we started from the CF structure for the J = 15/2 
ground-state multiplet initially derived by Ishikawa et al. 
for [DyPc2]~-TBA + ii2. Then we have varied the splitting 
among the levels until we found the best fit of the exper- 
imental data. Two sets of possible solutions were found 
to fit reasonably well %(T)T data. However, only one of 
them yielded a splitting between the lowest lying energy 
levels around 65 K, the value derived from AC suscep- 
tibility measurements (see later on). The corresponding 



Pm,m±i — C- 



mil 



E„ 



: {E m±1 -E m )/T _ I 



(5) 



Now the temperature dependent muon relaxation rate 
can be correctly fit using EqsOJ[5] (Fig. [TU]) by consid- 
ering the CF level splitting previously estimated from 
XsT fit and by using three different spin-phonon con- 
stants Ci and C 2 0, while C 3 ~ 2000 Hz/K 3 , which 
should be associated with the transitions \ m = ±11/2) <+ 
\m = ±13/2), \m = ±9/2) H |m = ±11/2) and 
jm = ±15/2) f> |m = ±13/2), respectively. We re- 
mark that C increases with the energy jump involved in 
the transition. This suggests that the relevant processes 
driving the transitions preferentially involve high energy 
vibrational modes. 
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FIG. 10: Temperature dependence of the muon longitudinal 
relaxation rate in [DyPc2]° for H = 1000 Gauss (circles). The 
line is the best fit according to Eqs. ((31 . 



In Fig. fTU)) one can notice that the best fit according 
to Eq|3] reproduces the muon relaxation data very well 
for T > 25 K. Since in this model only the thermally 
activated spin excitations are concerned, the behaviour 
of A(T) below 25 K should be probably ascribed to dif- 
ferent processes. For instance also in xsT a discrepancy 
from the theoretical calculation is observed at very low 
T. In fact, for T < 4 K xsT abruptly upturns, possi- 
bly signaling the onset of intermolecular correlations. In 
order to clarify this point, /iSR experiments have been 
performed down to very low T in [DyPc2]° and [TbPc2]° 
compounds, as shown in Fig. 6. As noted earlier, while 
the muon relaxation rate is nearly T-independent below 
1 K for both complexes (at least down to 50 mK), the 
very weak maxima located at T = 0.2 K and 0.7 K for 
the Dy and Tb samples, respectively, may in fact signal 
an onset of these intermolecular correlations. 

Now we turn to the comparison of the /iSR and AC 
susceptibility results. The maxima in x" Ixs shifting to 
lower temperature upon decreasing the irradiation fre- 
quency clearly indicate a progressive slowing down of the 
dynamics on cooling. In the case of a monodispersive dy- 
namical relaxation mode, at the peak temperature (T m ) 
the bulk magnetization relaxation time r c (T m ) matches 
the inverse of the angular frequency u) of the applied os- 
cillating field, according to the expression 



Xswt c 



1 



(6) 



We have checked the validity of this expression by per- 
forming frequency scans at a fixed temperature and found 
that in the explored T range the above expression was 
satisfied. The linear relation of in(r c (T m ) • ) to 1/T 
(Fig. [TT]) indicates that the Orbach process is domi- 
nant in the high temperature range and that the Ar- 
rhenius law t c = Toexp(A/T) for the correlation time is 



obeyed. From the fit of [TbPc 2 ]° data in FigfTTI one esti- 
mates a value for the energy barrier A ~ 750i-T, slightly 
smaller than the one deduced from muon relaxation rate 
measurements^ in the same sample. On the other hand 
it should be remarked that this value is much larger than 
the one reported by Ishikawa et al. (A ~ 590 K) in the 
same nominal compound^ The fit of [DyPc2]° data in 
Fig llll gives a much smaller activation energy A ~ 65 K, 
which corresponds quite well to the one obtained from 
the analysis of the static susceptibility data in Fig.([T|), 
and provides an estimate for the separation among the 
lowest energy levels of [DyPc2]°. This value is much dif- 
ferent with the barrier estimated from /iSR relaxation 
at higher temperature, which yields an estimate for the 
energy separation among the high energy levelsJ^ 
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FIG. 11: Natural logarithm of the spin correlation rate 
against the inverse of the x" /xs peak temperature for 
[TbPc 2 ]° (H = 1000 Gauss) (top) and for [DyPc 2 ]° (H = 9000 
Gauss) (bottom). The red lines correspond to the best fits for 
an activated behaviour for a barrier A = 750 K for [TbPc2]° 
and A = 65 K for [DyPc 2 ]°. 

In order to better compare the T-dependence of the 
correlation time derived by means of AC susceptibility 
and /iSR in [TbPc2]° we have derived, via Eq(51 the 
behaviour of r c estimated from AC susceptibility mea- 
surements and compared to the data reported in Fig. 
([7]), derived with /zSR. As it is shown in Fig. (fT2"]) the 
T-dependence of r c derived by the two techniques at dif- 
ferent magnetic fields overlap rather well above 45 K. At 
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lower temperatures a plateau is evidenced by both tech- 
niques and associated with tunneling processes among 
the |m = ±6) low-energy levels. It is noticed that the 
plateau lies at different values depending on the mag- 
nitude of the applied field and on the technique. As it 
has been pointed out in Refll2l the magnetic field leads 
to a Zeeman splitting of the two-fold degenerate ground- 
state and causes a reduction of the tunneling rate. Al- 
though this justifies the progressive increase of t c with 
magnetic field observed by each technique, it cannot ex- 
plain the significant difference in the r c values deduced 
by /zSR and AC susceptibility in the low-temperature 
tunneling regime. Thus, while both techniques probe 
the same high temperature activated dynamics driven by 
spin-phonon coupling, the low-temperature tunneling dy- 
namics probed at the microscopic or at the macroscopic 
level are different. This suggests that there is some cor- 
relation among the magnetic moments in the different 
molecular units which cause fluctuations which do not 
affect the total magnetization. This coupling should lead 
to flip-flop like fluctuations^ yielding a magnetic mo- 
ment flip from \m = +6} to \m = —6} on one molecule 
and the opposite flip in the adjacent molecule. These 
processes, which cannot be accounted for by a direct 
dipolar coupling among Tb 3+ spins since they produce 
a change Am = ±12, do not yield a net variation in the 
total magnetization and, hence, do not contribute to the 
AC susceptibility. In other words, at low-T /iSR probes 
T2-like processes involving the fluctuations of Tb 3+ mo- 
ments while AC susceptibility probes only those Ti-like 
processes yielding a net variation of the macroscopic mag- 
netization. The precise nature of these fluctuations, sim- 
ilar to the ones involving phonon trapping in NilO,— still 
has to be clarified. 
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FIG. 12: Temperature dependence of r c , deduced from Eq.© 
for a static field of 0, 500, 1000, 4000 Gauss in [TbPc 2 ]° com- 
pound. Solid circles correspond to t c temperature dependence 
as derived from fiSK analysis. 



Finally, we turn to the discussion of transport prop- 



erties. As it is shown in Fig.©, the temperature de- 
pendence of the resistivity is characterized by a rather 
small energy barrier, around 11 meV. Now, in [TbPc2]° 
molecular orbital calculations^ show that there is one 
unpaired electron in the ai highest occupied molecular 
orbital (HOMO), basically involving orbitals from the 
carbon atoms at the center of the Pc rings. If there is a 
sufficient overlap between the ai orbitals of the adjacent 
molecules, which are stacked forming chains, electron de- 
localization and a metallic behaviour should be attained. 
However, in those Pc-based systems the electron corre- 
lations are significant, particularly at half-band filling 22 , 
and the on-site Coulomb repulsion U can overcome the 
hopping integral t among adjacent molecules leading to 
a Hubbard insulating behaviour. Hence, the activated 
behaviour of the resistivity should be ascribed to corre- 
lation effects rather than to the band structure, which 
should typically lead to much larger activation energies. 
The band structure and, accordingly, the hopping inte- 
gral t appear to sizeably depend on the local structure, 
in particular, on the buckling of the Pc ring plane and 
on the rotation of adjacent Pc rings, hence it is possible 
that at high T the crossover from a negative to a posi- 
tive dp/dT with increasing T, might be ascribed to those 
structural effects. 

In conclusion, we have clearly shown that both AC 
susceptibility and /iSR probe a high T activated spin dy- 
namics in neutral [LnPc2] molecular magnets, character- 
ized by quantitatively similar correlation times. At low 
temperature both techniques exhibit a flattening of the 
correlation time which is associated with tunneling pro- 
cesses. Nevertheless, from /zSR one estimates correlation 
times which are one order of magnitude shorter than the 
ones derived by AC susceptibility measurements. This 
discrepancy between microscopic and macroscopic tech- 
niques indicates that there must be some correlation 
among Tb 3+ moments, which cause fluctuations which 
do not contribute to the uniform macroscopic suscepti- 
bility. The analysis of the static uniform susceptibility 
and of the /iSR relaxation rates in [DyPc2]° have allowed 
for a better characterization of the CF splitting of the 
J = 15/2 multiplet and of the spin-phonon couplings 
driving the activated dynamics in this compound than 
was reported in Refll3l. Finally, we have reported the 
behaviour of the resistivity in [TbPc2]°, where a low T ac- 
tivated behaviour, possibly arising from electron-electron 
Hubbard-like correlations, is evidenced. 
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